Deformation of geometry 
and bifurcations of vortex rings 



James Montaldi* & Tadashi Tokieda^ 
October 23, 2012 

Abstract 

We construct a smooth family of Hamiltonian systems, together with a family of 
group symmetries and momentum maps, for the dynamics of point vortices on surfaces 
parametrized by the curvature of the surface. Equivariant bifurcations in this family 
are characterized, whence the stability of the Thomson heptagon is deduced without 
recourse to the Birkhoff normal form, which has hitherto been a necessary tool. 

Introduction 

The present paper introduces one geometric idea and implements it in one dynamical prob- 
lem. 

Here is the problem. On the Euclidean plane, a ring of identical point vortices shaped 
as a regular «-gon is a relative equilibrium, in that it spins while keeping the same shape. 
Is this solution stable if we perturb the initial shape away from the regular n-gon? When 
n < 7, linear stability analysis (first carried out by Thomson [22]) concludes that the so- 
lution is stable; when n > 7, it is likewise unstable. But when n = 7, degeneracy makes 
linear analysis inapplicable and prevents us from concluding. Is the regular 7-gon stable 
or unstable? This is known as the Thomson heptagon problem. It has been answered in 
the affirmative, see [16, 11, 20], although as pointed out in [11] the argument in [16] is 
incomplete. Indeed part of our approach could be viewed as completing the argument of 
[16]. 

The spirit of the approach goes back to Poincare. In the theory of dynamical systems, 
the simplest solution methods to problems require some nondegeneracy condition (nonzero 
determinant, nonresonant frequencies, . . .). When the problem, call it V, is degenerate, we 
have to mobilize heavy machinery (cf. [20] for a proof with recourse to the Birkhoff normal 
form that the Thomson heptagon is nonlinearly stable, as well as for historical details). But 
there is an alternative approach. Embed V = Vo in a parametric family of problems and 
deform it away from degeneracy. The problem V can become tractable when regarded as 
lim^^^QVi, if we happen to understand well enough the bifurcations that occur in such a 
deformation. Thus, this approach trades one machinery for another, of bifurcation theory. 
The point is that the latter sometimes sheds an unusual light on the problem compared to 
the former. 

In the classic applications of this idea, people deform the dynamical system by adding a 
perturbation term in X. This, however, we cannot do in our special instance of the heptagon 
problem if we wish to preserve the hydrodynamic motivation: as long as we are studying 
the dynamics on the plane, it makes Uttle physical sense to tamper with the Hamiltonian. 
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We therefore deform not so much the dynamical system but rather the phase space on 
which the system evolves. A deformed choice of the phase space fixes canonically, by the 
hydrodynamic motivation, a deformed Hamiltonian formalism. Explicitly, we take X to be 
the Gaussian curvature' and deform the original plane to A > (family of spheres) and to 
A < (family of hyperbolic planes). Corresponding to this family of surfaces parametrized 
by A, we must write a whole parametric family of Hamiltonian systems for point vortices: 
a family of symplectic (Kahler) forms depending on A, a family of symmetry groups and 
momentum maps depending on A, a family of invariant Hamiltonians depending on X — 
all dependences arranged to be smooth. We do this in section 1. In section 2 we carry 
out the stability analysis for the parametric family. Bifurcations are characterized, and the 
nonlinear stability of the heptagon is deduced, in section 3. 

The idea of deforming the geometry underlying the dynamics of point vortices, in par- 
ticular as a route to a better understanding of the Thomson heptagon problem, arose during 
an evening conversation between the two authors in Peyresq, in the summer of 2003. We 
have since discussed it in seminars and conferences, and part of it has leaked into the lit- 
erature [1]. We set down the full story in the present paper The stability of the Thomson 
heptagon is stated below as Corollary 7. 

1 Smooth family of geometries 

In this section we describe the family of surfaces of constant curvature, containing the 
hyperbolic planes {X < 0), the Euclidean plane (X = 0), and the spheres {X > 0). Each is 
a homogenous space, i.e. an orbit of its group of symmetries, and for the different signs of 
X we describe the different groups. We begin with the 1 -parameter family of Lie algebras 
0;l (rather than groups) in section 1.1, and in order to obtain the corresponding family of 
surfaces A^;^, we shift the usual linear action of these Lie algebras to obtain affine hnear 
actions. In 1.2 the geometry of Ai^ is described; among other things it has curvature 4X. 
The Hamiltonian for point vortices on A^;^ is based on Green's function for the Laplacian 
on A4^, and we meet three possible choices of Green's function according to choices of 
the 'boundary condition'. In 1.3 we comment on the implications of the different choices 
of Green's function. 

1.1 Lie algebras 

Let A e M. On M? with coordinates {x,y,u), consider the family of metrics 

ds^ ^dx^ + dy^ + Xdu^. (1.1) 

The Lie group of linear transformations preserving this metric will be denoted SO{qi), 
where qi = diag[l, 1, A] is the metric tensor For the Lie algebra, we have X E 5o{qi) if 
and only if X^qi +qiX = 0. A basis for so{qi) is 

/O \ / A\ /O -1 0\ 

Xi= -A , ^2 = , ^3= 1 . (1.2) 

yo 1 / \-i 0/ yo 0/ 

(Strictly speaking, the third column of X, is arbitrary when A = 0, meaning the family of 
all automorphisms of (1.1) is not flat. We are picking a component which is a flat family 
over A.) This basis satisfies the commutation relations 

[^1,^2] — ^^3, [-'''2,^3] =^^1, [X3iXi]~X2. 

From now on, we shah abbreviate Gi = SO{qi) and — so{q;^). 

' Actually formulaic convenience leads us to take 4A to be the Gaussian curvature. 
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The Lie algebra is isomorphic to so(3) for A > 0, to se(2) for X =0, and to s[(2) for 
X <Q. Indeed, for A 7^ 0, the standard commutation relations are recovered by rescaling 
the basis to {|A|-i/2Xi, \X\-^/^X2,X3}. So 

r S0(3) ifA>0 
^ < SE(2) if A = 
[SL{2,M) ifA<0. 

It is seen from the commutation relations that in the adjoint representation of Qi the 
basis elements are represented by adx, = —Xj ; in other words, if Y,] ojXj G g;^ is written as 
a vector u = (ai 02 ^3)^, then adx^ (u) — —Xju. In the coadjoint representation the basis 
elements Xj <E Qx are represented by the matrices Xj themselves. Thus the original from 
which we started may be naturally identified with q*^ . 

Afflne action Whereas the coadjoint action defined by the matrices Xj depends contin- 
uously on A, it is not possible to track a single orbit continuously as A crosses 0. Yet, in 
what follows, we wish to do just that. To this end we must shift the linear coadjoint action 
by a translation, making it an affine action. 

The affine action of the Lie algebra is given by 

X-pi^Xpi + x{X), (L3) 

where X/i is the linear action part (matrix times vector) and 

(-b/2 

T(flXi +bX2 + cX^) = fl/2 

V 

is the translation. The orbit we track is the one through the origin, cf. section 1.2. 

Remark 1. Our translation t, a function of the element of g;^, is a 1-cocycle taking values 
in Q*^, a symplectic cocycle in Souriau's terminology [21] because the matrix of T is skew- 
symmetric. It is known that every cocycle is exact when the group is semi-simple, and our 
Gi is semi-simple for A 7^ 0. Here T = 5((0, 0, 1 /2XY), since by definition 5{pL){X) = 
— coadxjU = —Xjj.. The natural invariant Poisson structure on = with the cocycle T 
is given by (cf. [15, 21]) 

{/, g}{fi) - {fi, [df{fi),dg{m - (T(d/(M)),dg(Ai)) (1.4) 

under the identification d/(/x),dg(/x) G (fl^)* — 0a,- The Casimir for this Poisson structure 
is +y^ + Xu^ — u, so level sets of this function are the orbits of the shifted coadjoint action 
(1.3), of which we shall take advantage below. For the record, the Kostant-Kirillov-Souriau 
symplectic form on the affine coadjoint orbits is given by the same formula, 

n^(u,v)-(M,[^,J?]>-(T(<^),rj), 

where u = coad^/i and v = coadjj/i. 

1.2 Surfaces 

Now consider the family of quadratic surfaces through the origin in R^, 

x^+y^ + Xu^~u = 0. (1.5) 

When A > 0, this looks like an ellipsoid with centre at {x,y,u) = (0,0, 1/2A). With the 
metric (1.1), however, this ellipsoid-looking surface is in fact a sphere of radius l/2\/A. 
Its Gaussian curvature is 4A. 
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X > (sphere) 



A < (hyperboloid) 



Figure 1.1: Geometries in the 1 -parameter family 



When A = 0, (1.5) defines the paraboloid u=x^ +y^, and the metric is the usual metric 
on the jcy-plane lifted to the paraboloid by orthogonal projection, so is of curvature 0. 

When A < 0, the metric (1.1) becomes Lorentzian, but restricted to either sheet of the 
2-sheeted hyperboloid defined by ( 1 .5) it induces the hyperbolic metric of constant negative 
curvature 4A; we consider just the 'upper sheet' that passes through the origin, see Figure 
1.1. 

We refer to the surface (1.5) with metric induced from ( 1 . 1) as A^;^ . It is easy to check 
that Ail is invariant under the infinitesimal action (1.3), and is therefore an orbit of the 
affine coadjoint G;^ -action on g^. 

To create a uniform coordinate system on A^;^, we use stereographic projection on the 
jcy-plane, centered at the point (0,0, 1/A) (where M.x intersects the M-axis, besides the 
origin); for A = this is the orthogonal projection. We also identify the xy-plane with C, 
via z~x + \y. The map inverse to the projection has the formula 

The domain of this map is {z G C | 1 +A|zp > 0}, which is the entire plane if A ^ 
and a bounded disc (Poincare disc) if A < 0. For the sphere, the equator corresponds to 
|zp = 1 /A, while the point antipodal to z is — 1 / A|. 

The metric on the surface M.x induced from that in (1.1), in terms of the complex 
variable z, is 

d/ = -i|dz|2, 

where a = 1 + A |zp, a notation we shall use throughout. The circle |z| = r in C maps to a 
circle of radius a on A^;^, where^ 

{ -^taxc^ [rVl^ ifA>0 
r ifA=0 (1-7) 

tanh- 1 {r^/^) if A < 0. 

Pulling back the vector fields Xj of the Lie algebra (or rather their affine variants shifted 
by t) via the stereographic projection yields 

^,{x,y)^\{2Xxy., 1-A(x2-/)), ^2{x.y) ^ ~\{\ + l{x^ -y^), iXxy) , 

'Despite three formulae, a is a single analytic function of nA, with series expansion a = r— l^r'X + ^f^X^ — 
^r'A^ H convergent for \r-X\ < 1. 



or in complex variables 



or in polar coordinates 



^i+i^2 = ^ei«^^5e + ^ie'V5,, ^i = de. 

Symplectic structures Up to a scalar multiple, there exists a unique SO(^;L)-invariant 
symplectic form on We choose the scalar so that 

2 i 

Q.l = ^(ixAdy= ^dzAdz. (1.8) 

The choice of scaling is such that the sphere Aii of radius l/lVX acquires symplectic 
area 2nX^K With respect to the basis {Xi,X2,Xi} for the Lie algebra, the momentum map 
takes the form 

h{z)^^izM^). (1.9) 

This coincides with the inclusion Aix^M? given in (1.6), which shows that Q.x coincides 
with the KKS symplectic form on the affine coadjoint orbit. 

Green's functions The metric (1.1) on induces the metric on Aix- In terms of the 
uniform coordinate system (1.6), the metric tensor is (7^^diag[l, 1]. The Laplace-Beltrami 
operator on is 

The 2-point Green's function for this operator is 

G(z;w) = log |z-w|2. (1.10) 

This satisfies A-G — for z^w and has a logarithmic singularity at z = w. When we regard 
the plane as a model for (most of) the sphere, G has another singularity at z = oo. 
An alternative Green's function is 

This satisfies A,Gi (z; w) —0 for z ^ {w, — I /Xw} and has a logarithmic singularity at those 
excluded points (which are antipodal to each other), but is regular at z = oo. 

The Green's function that is usually used on the sphere is the 'log of Euclidean dis- 
tance', whose expression after stereographic projection is 

|Z- w|2 

^^^'"^='°g (l+A|z|2)(l+A|wP) - ^'-'^^ 

Away from the pole at z = w, this satisfies A,G(z;w) = — 4A, so it is not, if we go by the 
book, a Green's function. This function is regular at z = °°, and also has a well-defined 
limit as w — > oo, namely G(z;°°) — — log (7 (up to an additive constant of —log A), which 
checks against Aj.(— log a) — — 4A. 

In the next section we comment on the differences among these Green's functions in 
the context of the dynamics of point vortices. 
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1.3 Hamiltonians for point vortices 



We recall how the Hamiltonian formalism for the the dynamics of point vortices works. 
Let (m, v)^ be the velocity field of an inviscid, incompressible flow on a domain D C 
~ C. The incompressibility + = implies the existence of a stream function 

V/ : D ^ R such fliat 

u=^\j/, v=-^\i/. (1.13) 
oy ax 

The curl of the velocity ' is — — —Axj/. The boundary condition for an inviscid flow 
is that (m, v)^ be tangent to dD everywhere, equivalently that every connected component 
of dD be a level set of y/. The total circulation along all the boundary components is, by 
Stokes's theorem, 

/ udx + vdy^ / -Axjfdxdy, (1-14) 
JdD Jd 

the total curl present on D. 

Now consider a model situation where the flow is generated by a curl concentrated at a 
singularity zq = -^o + iyo and the circulation around that singularity is Itzk: 

-Ava(z) = 2kk5{z-zo). 

We recognize that, up to the sign and a scalar coefficient, \j/ is Green's function. We say that 
we have a point vortex of vorticity K at zo- In the situation where we have point vortices 
of vorticities Ki,...,k„ at zi , . . . ,z„, each vortex moves carried by the sum of the flows 
generated by all the other vortices. The equations (1.13) show that the dynamics then is 
Hamiltonian. 

The theory is written analogously on any domain of any Riemann surface. In particular, 
on the surfaces Aii discussed above, the Hamiltonian is 

(Zl , . . . , ZAt) = - £ K-i K-jG;t (z,-; Z/) . 

><j 

Note an unusual feature of this Hamiltonian system: unlike in classical mechanics, the 
phase space here is not the cotangent bundle of anything, but rather the n-fold product 
M.X X • • • X A^A (minus the diagonals if we want a priori to avoid collisions) with the 
weighted-sum symplectic form K\Q.x • • • Kn^x . 

For A > 0, A^;l is compact without boundary. In this case, the fact of nature (1.14) 
forces the total vorticity to be zero: Y.] = 0. This, in principle, bans placing a lone point 
vortex on A4x for A > or on any closed Riemann surface. The dodge around this ban, 
favoured in the Uterature, is to impose a constant background vorticity 

—sum of circulations / area — —X ^ Kj, 

J 

which results in (1.12). For a (geo)physical example, a rigidly rotating sphere entails such 
a background vorticity. But even if background vorticity dodges around A > 0, continuing 
it to A ^ gets us into trouble, for on these noncompact surfaces, background vorticity 
imparts an infinite amount of energy to the flow. So for A ^ the family (1.10) seems 
preferable. 

However, welding together (1.12) for A ^ 0, (1.10) for A < has a decisive defect: the 
resulting family is not smooth in A. There are three options for having a smooth family. 

1. Use (1.10). This costs postulating an immobile vortex of vorticity — the North 

Pole for A > 0. 

^The minus sign makes —A a positive operator. But we sliall be casual about the sign and use +A as well as 
—A. All that the casualness causes is to reverse the direction of the flow. 
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2. Use (1.11). This costs introducing 'counter-vortices' at antipodes to the zjs for A > 

(recall that the point antipodal to z is — l/Xz). 

3. Use (1.12). This costs infinite energy for A ^ 0. 

In the planar case A = all three Green's functions (1.10), (1.11), (1.12) agree. 

In this paper, we opt for the family defined in ( 1 . 1 2) with the constant background vor- 
ticity, because the infinite energy of a tame flow would not shock any fluid dynamicist — 
flows on the plane uniform at infinity and such are handled routinely — whereas the smooth- 
ness of the family, without postulating extraneous objects, is essential for us. In section 4 
we sketch how the analysis can be adapted to the other two options. 



2 Nondegenerate analysis of vortex rings 

In this section we study a ring of n vortices with identical vorticities K, shaped as an n-gon 
on A^;^. We evaluate the Hamiltonian, the momentum map, and the augmented Hamilto- 
nian at the regular ring in section 2. 1 . We obtain the Hessian of the augmented Hamiltonian 
and its spectral data in 2.2. In 2.3 we construct the symplectic slice and say what we can, 
as far as this linear analysis goes, about the stability of the ring. The answers will depend 
not only on n but also on A . 

Recall that the augmented Hamiltonian is given hy H — coJ, where J ^ Jx is the con- 
served quantity coming from the rotational symmetry; it represents the Hamiltonian in a 
frame rotating with angular velocity o. Its critical points are therefore equilibria in the 
rotating frame, i.e. relative equilibria, and the Hessian restricted to the symplectic slice 
determines the stability. 

2.1 Regular ring 

Let 7 be a primitive n th root of unity. When the vortices are placed at the vertices Zj = 
a regular n-gon of radius r, the Hamiltonian Hx takes the value 

where Gx — logop;^ with Pxir^rY^) = |1 — 7-'pr^/(l + Ar^)^. In view of the identity 



.1 



hence 



, , 7, n(n- 1)k-2, P- 



(1+Ar2)2 

recall the momentum map given in (1.9). The first component vanishes for these regular 
rings, while the value of the second component at the regular ring is 



(Notice that the full momemtnum map is typeset in bold, while this component is not.) 
Ignoring the constant, the augmented Hamiltonian then takes the value 

«(«— 1)K'2 ^2 ^^^2 



Stt ^(1+Ar2)2 "'1+Ar2' 
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which admits a critical point at r = ro 7^ if and only if 

(n-l)K 

(0 = (0o^-^^ (2.1) 

This is the angular velocity of the regular ring. (It is a Uttle surprising that (Oq is even in A.) 

2.2 Hessians 

We continue with the system of n identical point vortices, all with vorticity k. In the 
uniform coordinate system on Mx, the Hamiltonian with Green's function (1.12) is 

where = 1 +A|zj|2 (j = !,...,«). The augmented Hamiltonian then is 
Hf^izi, . . . ,z„) ^ Hi{zi, . . . ,Zn) - coY, 1^-^- 

HI ^/ 

We saw above that this is critical at Zj = roe^"^^/" and CO = COo as in (2.1). The entries in 
the Hessian are 



52 ^ _ k2 2 ^ ^ _ )f2 " 

52 



At the critical point, and with O) = COq, we have 



J^^k 4;r(l-cos; 
//a =0 (V;,^). 



(2.2) 



where 

(so that C7 + a = 2). We used the identity, valid for ^ £ ^ n, cf. [8]: 

-^os(2^^ 1 

l-cos(2;r;/n) 6^ ^ ^ ^ 

Eigenvalues of the Hessian The Hessian matrix d^H^ is block-diagonal, with two nxn 
blocks both of which are symmetric circulant, so the eigenvalues can be written down at 
once. 

Following the notation in [ 13], for £ — 0,1, ... ,[11 /2J define the Fourier tangent vectors 

C/^'' = aj^^^+ipP = Yexp{-2m£j/n)5rj 

n (2-5) 
C^'' = a(')+i^f = i^i:exp(-2;ri £;■/«) 50,. 
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Here 5rj denotes the unit tangent vector in the -direction, and similarly for 50j/ro. For 

£ = 0, n/2, we have jS/^'' = jSg^'' = and the ^'^'s are real. The as and j3s form a set of 2n 
linearly independent vectors, forming a basis for the tangent space. For the record, 

5rj^-j^exp{2mij/n)i;l'\ 

and similarly for 5 9j . The span Vi = af^ , a^^^ , jS/^' , j3g^'' ^ is a subspace of Fourier modes; 
for £ = 0,n/2 they are 2-dimensional, while for all other indices £ they are 4-dimensional. 

As each block of d^Hi is circulant as well as symmetric, i^/^' and ^g^' (or rather their 

(f) (t) 

real and imaginary parts) are the eigenvectors. The eigenvalues £r and (of course 
real) are found to be 

= A + ^((„2-l)_6^(«-^)) 

i2A 



(2{n-l)'-±^~£in-£)^ (2.6) 



4'' = 

where A is given in (2.3). In case £ = 1, these simplify to 



^2 



A-nrQG'^ An 

which are both strictly positive. 

2.3 Symplectic slice 

Not all the eigenvalues are relevant to stability. First, those that are zero because they corre- 
spond to directions along the group orbit should be discarded. Second, those corresponding 
to directions transverse to the level set of the conserved quantities should be discarded, too. 
This is the process of reduction: restrict to KerdJ (where i — ix is given in (1.9)), then 
take the complement to the group orbit in that kernel. The resulting space is called the 
symplectic slice, which we denote by M. 

We find KerdJ using the Fourier basis above. Identify with C x R. Then 

while dj vanishes on all other Fourier tangent vectors. Write V[ = Vi n KerdJ. Then 

aoi'^-afi'i\ (7j3/'' + 5a<^'). (2.8) 

KerdJ is spanned by the V^s with £> \,V[, and i^g"'. The subspace generated by tl}p is 
tangent to the group orbit (an infinitesimal rotation about the origin), so must be discarded. 
Finally the symplectic shce is 

L«/2J 

U = y[®^Vi. (2.9) 

The relevant eigenvalues are therefore , fig for £ ^ I. Now, with respect to the basis 
for Vj given in (2.8), the restriction of the Hessian to V[ is a scalar multiple of the identity. 
So it has a double eigenvalue 



^ I a\ (2.10) 
2 2rQ Altr^ 
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which is strictly positive unless Xr^ = 1. However, Xr^ = I corresponds to the equator 
on the sphere, at which the momentum value is left fixed by all of S0(3). In this case J\f 
drops in dimension, V[ no longer lies in J\f, and the corresponding eigenvalue e[ becomes 
irrelevant to the stabiUty of the relative equilibrium. 



Stability The relative equilibria in question are rotating rings; thus they are periodic tra- 
jectories. The precise sense of stability we are adopting is like the ordinary one of Lyapunov 
stability, but in terms of -invariant open sets, where jj. is the momentum value at the tra- 
jectory and Gjx is the stabilizer of jj.: in detail, for every G|j -invariant neighbourhood V of 
the trajectory, there exists a G^-invariant neighbourhood U CV such that every trajectory 
starting in U remains in V for all time. 

We can make do with a coarse criterion: if the restriction of the Hessian of the aug- 
mented Hamiltonian to the symplectic slice is positive-definite, then the relative equilib- 
rium is stable in our sense. A finer criterion is: if this Hessian is merely non-negative but 
the augmented Hamiltonian admits a local extremum at the relative equilibrium, then the 
relative equilibrium is still stable [17, Theorem 1.2]. 

Among the relevant eigenvalues, e^^' > for £ ^ 1 and ej > 0. It remains to check 

the sign of for £ ^2. From the expression in (2.6), it is clear that the least eigenvalue 
occurs for £— [«/2j . This is the criterion we have been after: the relevant eigenvalues are 
all strictly positive if and only if 



1+AM 1 



'0 > 



(1+Ar2)2 2(«~1) 



(2.11) 



The left-hand side is unbounded as a function of ro if A < 0. The right-hand side is a 
strictly increasing, unbounded function of « for « ^ 3. Hence, on one hand, for fixed A, ro 
a value of n exists beyond which the inequality fails; on the other hand, if A < 0, for fixed 
n, by enlarging sufficiently close to — 1 /A (its supremum on the hyperbolic plane), we 
can ensure the inequality holds. We conclude with a result which for spheres is due to [2] 
(see also [13]): 



Tlieorem 2. On A4i, the relative equilibrium ofn identical point vortices in a regular ring 
of radius a is Lyapunov- stable if (2.11) is satisfied, where a and ro are related by ( 1.7). In 

'^0 



particular, in the hyperbolic scenario A < this ring is stable if Xr^ is sufficiently close to 



-1. 

Remark 3. Physical intuition confirms that a spherical surface destabilizes a ring whereas 
a hyperbolic surface stabilizes it. Think of a vortex of the ring. As A gets positive, the 
adjacent vortices remain relatively near while the diametrically opposite vortices become 
relatively far, so the former, which tend to knock our vortex perpendicularly to the ring, 
exert more influence than the latter, which tend to slide our vortex tangentially to the ring. 
As A gets negative, the effects are felt the other way round. 

Remark 4. If the inequality in (2. 1 1) is reversed, then the ring is linearly unstable. This is 
because the V[„/2j mode will have real eigenvalues. Indeed, for n even y„/2 is of dimension 
2 and the Hessian is indefinite, which suffice to conclude that the eigenvalues are real. For 
n odd V(„_i)/2 is 4-dimensional, so having an indefinite Hamiltonian does not imply the 

eigenvalues are real. However, the negative eigenspace is spanned by ai^' , jS/^^ which is 
Lagrangian, and then the realness of the eigenvalues follows. 

If the parameters are such that (2. 1 1) is an equality, then the stability is not determined 
by linear analysis. This determination is what our deformation plus bifurcation approach 
achieves, in sections 3.5.1 and 3.5.2. 



We now spell out the conclusions concretely for each of the three geometries. See 
Figure 2. 1 and the table of bifurcation points below. On the spheres A > there are two 

.2 



bifurcation points: the value of Ar^ listed in the table below and its reciprocal 
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-1 1 loga 

Figure 2.1: Ranges of stability (solid lines) and instability (dotted lines) for the dimension- 
less quantity log a = log(l + Atq) 



n 


4 


5 


6 


7 8 


9 


10 


11 


12 


13 


H 


0.268 


0.172 


0.0557 


-0.0627 


-0.101 


-0.143 


-0.172 


-0.202 


-0.225 



Plane There is an obvious scale-invariance, and the stability/instability of the ring is in- 
dependent of the radius a. We recover J.J. Thomson's original result [22] that the ring is 
stable when n < 7 and unstable when n > 7. When n = 7, we have the so-called Thomson 
heptagon, whose stability is not determined by linear analysis (but cf. Corollary 7). 

Spheres When the ring is always unstable. When n <1, there is a range of a over 

which it is stable. Example: when n — 6, it is stable for (Atq — 2)^ > 3 which translates, 
bearing in mind A > 0, into 

2 2-V3 2 2 + V3 

ro < — — or Kq > — — . 

The two inequalities correspond to neighbourhoods of the South and North Poles respec- 
tively. For n ~3 the ring is always stable, at any radius a. The linear stability of rings of 
vortices on the sphere was first studied by Polvani and Dritschel [18], and the full nonlinear 
stability in [2] — see also [13] for more details. 

Hyperbolic planes When « ^ 7, the ring is always stable. When « > 7, it is stable (for a 
given A) for a sufficiently large. Example: when n = 15, it is stable for 

2 2-V3 

3 Bifurcations across the degeneracy 

To understand the bifurcations in detail, it is imperative to exploit the symmetries of the 
system, which for the ring of point vortices is the dihedral group. We begin section 3.1 
by setting up the dihedral symmetry of the system, and then in 3.2 list the bifurcations we 
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expect in Hamiltonian systems with dihedral symmetry. Following that, in 3.3, we state the 
main theorem (Theorem 6) on which of these bifurcations occur in the dynamics of point 
vortices; it leads to the nonlinear stability of the Thomson heptagon (Corollary 7). In 3.4 
we visit the geometry of the bifurcating rings, which enjoy less symmetry (spontaneous 
symmetry breaking) and are illustrated in Figure 3.4. Finally, 3.5 proves the main theorem: 
we perform the calculations needed to justify these 'expectations', and decide which of the 
expectations are actually realized. For the sake of completeness, section 3.6 summarizes 
results from [ 14] on bifurcations from the equator; they are not covered by the other results 
as the momentum value there is degenerate (in the sense that it is fixed by the entire group 
rather than just a 1 -parameter subgroup). 

3.1 Dihedral group action 

For the simplicity of language, we shall confuse M.^ and the uniform chart C of section 
1.2. Points and group actions on C should be interpreted as their lifts on A^;^ . 

The full system is invariant under the symmetry group Gx (depending on X) as in 
section 1.1. For every A, Gx contains rotations about the origin, and reflections in lines 
through the origin, together generating a subgroup of Gx isomorphic to 0(2). Consider 
now a system of n identical point vortices. A dihedral subgroup^ D„ C 0(2) x Sn acts on 
the phase space C" by 

c- (zi,...,z„) = (cz„,czi,...,cz„_i) ^^^^ 

where c — exp(27ri/n) is a cyclic rotation and m is a mirror reflection. If the points Zj = 
roexp(2;riy/n) are placed as a regular ring, then that configuration is fixed by this D„, and 
m acts as a reflection in the line passing through z„ . In case n is odd, all reflections in D„ are 
conjugate, whereas in case n is even, there are two distinct conjugacy classes of reflections, 
one consisting of those through opposite vertices of a regular n-gon (all conjugate to m), the 
other consisting of those through mid-points of opposite edges (all conjugate to m' = cm). 
This will come into play in deciding what bifurcating solutions appear. 

Since the D„-action fixes the ring, D„ acts on the tangent space to the phase space at 
this ring. The Fourier basis (2.5) for n point vortices is adapted to this action: 

and m' ■ C/^' = Q-27tie/n . ^ _^-2m(/n ^in^ Pqj. fjjg 2-dimensional subspace V/ 

of the symplectic slice JV (2.8), write 

C = (era/') - aj3^' + = aC/" +iaC^'). (3.3) 

Then (3.2) implies that c • C' = e-^""'/" m • C' = and m' ■ C = e-^'i'/" C'. 

As the parameter X varies, the eigenvalue e/^', £ ^ 2 in (2.6) may cross and may 
involve a bifurcation in the mode V(. In contrast e^'' and eQ^\ being strictly positive, never 
involve bifurcations and in particular the mode never bifurcates. 

3.2 Dihedral bifurcations 

The type of bifurcation expected in a symmetric Hamiltonian system is controlled by the 
group action on the generalized kernel of the linear system at the bifurcation point. This 
was first investigated by Golubitsky and Stewart [6], cf. also [4]. It follows from [6] that 

"^S,, is the group of permutations of tlie n vortices and D„ lias order 2n. 



12 



in a generic family of linear Hamiltonian systems, a pair of eigenvalues come together 
along the imaginary axis, collide at the origin, and split along the real axis. This splitting 
transition is indeed what happens in our problem, as seen from the expressions (2.6) for 
the eigenvalues. Part of the genericity hypothesis of [6] is that the generalized kernel be an 
irreducible symplectic representation, which is satisfied here. 

The greatest common divisor of n and i will be denoted by (n,^). We see from (3.2) 
that the cyclic subgroup Z(„ C D„ acts trivially on Vi. Consequently on Ve there is an 
effective action of D„/(„ ^j. It turns out that Vi is an irreducible symplectic representation of 
this group. Two cases are to be distinguished: £ = n/2(n even) when V„/2 is 2-dimensional 
with an action of D2 = Z2 x Z2, and £ ^ n/2 when the V^s for Q < £ < n/2 are all 4- 
dimensional. For bifurcations the modes 1 < £ alone are of interest to us. Much of 

the analysis of generic bifurcations with dihedral symmetry is found in [7], and although 
there they deal with general vector fields rather than with Hamiltonian ones, the conclusions 
turn out to be the same. Analysis of the gradient case is also in [3]. 

In the dynamics of point vortices, if the bifurcating mode is ^ 7^ [n /2J , then the Unear 
system"" at the relative equilibrium has real eigenvalues in the [«/2j mode, hence is unsta- 
ble. If £= [«/2j , then the bifurcation involves a loss of stability in the mode £. This means 
that in our analysis below, a local minimum corresponds to stable relative equilibria only if 
we are looking at the [«/2j mode. 

Bifurcations on Vc for £ = n/2 This is the 2-dimensional symplectic span 



with an action of D2 ~ Z2 x Z2. The kernel of the Hessian at the bifurcation point is the 
1 -dimensional subspace spanned by (}"^^\ with a Z2-action. Then the generic bifurcation 
is a Z2-pitchfork, which can be either sub- or super-critical; if subcritical, the bifurcating 
solutions are unstable, and if supercritical, they are stable (provided the original 'central' 
solution was stable). A normal form is given by the family 



'h.o.t.' stands for higher-order terms in x^,y^,M. The + sign in front of is justified by 
the eigenvalue in the ^g"'^^' -direction, which is always positive. The — sign in front of u 
is a choice, dictated by the fact that increasing kr\ makes a critical point pass from local 
minimum to saddle, as shown in Figure 2.1. The sign + or — in front of corresponds to 
supercritical or subcritical, respectively. See Figure 3.1(a,b) and the lecture notes [4] for a 
fuller discussion. In Remark 1 1 we explain why the bifurcations occurring here are in fact 
all supercritical, for all even n. 

Bifurcations on for £ 7^ n/2 Of these only £ — j{n — 1) involves stable relative equi- 
libria, other modes bifurcate only if the linear system already has real eigenvalues; that 
said, the other values of £ do involve the appearance of new unstable relative equilibria, so 
are of interest. Write k — n/{n, £). Then A; ^ 4 and Dj. acts effectively on Vf. The type 
of generic bifurcation we get depends on k. The 4-dimensional Vt is a direct sum of two 
Dyj-invariant 2-dimensional Lagrangian subspaces, on one of which the Hessian vanishes 
at the bifurcation point, on the other it is always positive-definite. A Di.-invariant function 
on such a space is a function of the fundamental invariants 




fu{x,y) 



Mx^±x'^+/ + h.o.t. 



(3.4) 



N{x,y)=x^ +y^, P{x,y) ^R&ix + xyf . 



'The vector field, not the Hessian. 



13 



(a) D2 supercritical pitchfork (b) D2 subcritical pitchfork 
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(c) D3 (transcritical) (d) D4 transcritical (|a| < |j3|) 




u 



(e) D<. pitchfork, k \ 
(with \a\ > 1)3 1 foryt = 4) 



Figure 3.1: Bifurcation diagrams for generic Dj. -bifurcations, 2. r — \/ +3^^, and u 
is the parameter as in the text. Solid lines refer to local minima, dashed lines to saddles and 
local maxima. For ^ ^ 3 each nontrivial branch corresponds to k solutions after applying 
the rotations from D^.. Figure (e) is drawn for a > 1, where a is the coefficient in (3.5); if 
0: < ~ 1, reflect the diagram in the r-axis. 



A generic 1 -parameter family of such functions is given by 

/„ = -MA'+aA'2_^i3P + h.o.t., (3.5) 

where 'h.o.t.' stands for higher-order terms in A^, P, u. Figure 3.2 shows the level sets of /„ 
for k — 3,. . . ,6, as u varies through 0. 

• If ^ = 3 or if A; = 4 with ji > a, then the bifurcation can be said to be transcritical, in 
that the bifurcating branches exist on both sides of the bifurcation point u — 0, and the k 
bifurcating points are all saddles (and hence unstable equilibria), each of which is fixed by 
a mirror reflection conjugate to m. See the bifurcation diagrams in Figure 3. l(c, d). 

• If A: = 4 and a > )3 or if A: > 4, then the bifurcation is like a pitchfork, in that all bifurcating 
equilibria coexist on the same side of the bifurcation point. But unlike the pitchfork, 2 types 
of bifurcating solutions appear, possibly with different stability properties; if k is even, then 
one has symmetry type (m) and the other (;«') . See the bifurcation diagram in Figure 3. 1(e). 

Remark 5. The finite determinacy and unfolding theorems of singularity theory guarantee 
that /() in (3.4) and (3.5) is finitely determined and h.o.t. may be ignored. If n ^ 3, then 
/o has codimension 1 provided j3 7^ 0, and /„ is a versal unfolding of /o, so that any de- 
formation is equivalent to it. If n ^ 4, then fo has codimension 2, and a versal unfolding 
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(a) D3 with « < 



(b) D3 with u > 




(c) D4 with \a\ < \p\ and 
u<Q 



(d) D4 with \a\ < \p\ and 
u>Q 



Figure 3.2: Contours of the generic 1-parameter family of D|(.-invariant functions (3.5), for 
k = 3,4 [produced with Maple, with a judicious choice of level sets] . The figures (a)-(d) 
are all transcritical bifurcations. Continued on next page. 



provided j3 ^ (and a ^ ±j3 when n = 4). The parameter v defines a topologically trivial 
deformation, i.e. v can be eliminated via a continuous change of coordinates rather than a 
smooth one; nevertheless this homeomorphism will be a diffeomorphism away from the 
origin, so critical points are preserved. The modulus v that arises when n = 4 is related to 
the cross-ratio of the 4 lines making up (0). 

3.3 Bifurcations of vortex rings 

We work with the parameter Atq > —1. Recall that ro is the radius of the vortex ring 
measured on C after the stereographic projection, and that it is related by ( 1 .7) to the radius 
a measured on Mi . The curvature case is A = 0, and in the spherical case A > the 
values Atq and l/Ar^ are equivalent as they represent antipodal rings on the sphere. We 
therefore let Xr^ vary in the range (—1, 1]. 

Now for the main theorem. Stability means Lyapunov stability modulo rotations (same 
as orbital stability in our situation). Instability means full spectral instability, i.e. at least 
one of the eigenvalues is real and positive. The ring of n = 3 vortices is always nonlinearly 
stable. 

Theorem 6. Let A. With Xr^ G (— 1, 1] oi a parameter, the regular ring ofn identical 
vortices undergoes the following bifurcations, illustrated in Figure 3.3: 

all n The ring is stable for Atq < b„, where b„ is the unique root in (—1,1] of^ 



IS 



fu,v 



uN+{a + v)N^ + ^P 




(3.6) 



6, 



•cf. (2.11). 
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(e) D4 with \a\ > \ f}\, u< (f) D4 with |a| > \p\, u > 






(i) Dfi with M < (j) De with w > 



Figure 3.2: (continued) Contours of the generic 1-parameter family of D|t-invariant func- 
tions (3.5), for k^4. These figures are all dihedral pitchfork bifurcations. 



n even As Xr^ crosses bn, the ring loses stability via a supercritical pitchfork bifurcation, 
and the (stable) bifurcating solution consists of a pair of n/2-gons with different 
values for the radius. 

n odd As Xr^ crosses b„, the ring loses stability via a supercritical bifurcation as depicted 
in Figure 3.1(e), to 2 types of relative equilibria, each with a line of symmetry. 

all n As XrQ increases further, the ring undergoes a sequence of bifurcations, one in each 
of the modes [n/2j > i> \; all the relative equilibria involved are unstable, and the 
bifurcating solutions have D(„ (ysymmetry.^ 

As a special case, we recover the following result of Kurakin and Yudovich [11] and 
Schmidt [20]. The calculation justifying it is the subject of section 3.5.2. 

Corollary 7. The Thomson heptagon is nonlinearly stable. 

Bifurcation values of XrQ The tables below spell out the values of Xr^ where the bifur- 
cations occur, for n ~ 6,7,8,9. They are found by solving e^^'' — (2.6) for Atq; the first 
values are those of b„ mentioned in Theorem 6. 

'We put Di = Z2 acting by reflection. 
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(2) 



(2) 



(3) 



(2) 





Figure 3.3: Bifurcation diagrams for the ring of n identical vortices for low values of n. The 
number in parentheses is the mode number bifurcating at that point. The black dot on the 
axis represents schematically the point where A = (the plane). X increases toward right. 
We do not know whether the bifurcating branches from the lower modes branch to the right 
or the left, though we believe they are as shown. The case n = 8, £ = 2 has an effective 
action of D4, so could be transcritical or pitchfork — we do not know which occurs. 





n = 6 






n = l 




mode 


1 = 3 


£ = 2 


mode 


1 = 3 


1 = 2 


Xrl 


0.056 


0.127 


Xrl 





0.101 





n = 8 




n = 9 




mode 


£ = 4 £ = 3 1 = 2 


mode 


e = 4- 1 = 3 


1 = 2 


Xrl 


"0.063 -0.033 0.084 


Xrl 


-0.101 -0.056 


0.072 



In all the tables, the bifurcation of the 1 = 2 mode occurs for A > (on the sphere); this is 
easily checked to be true for all n. 

3.4 Geometry of bifurcating rings 

At a bifurcation, the bifurcating mode controls the geometry/symmetry of the bifurcating 
solution. Points in all correspond to configurations with cyclic symmetry Z(„ C D„, 
which allows the D„-action on to factor through a D„/Z|-„^) ~ D^-action, where 

k = nl(n,l). 

Moreover, the bifurcating solutions are all fixed by a reflection in D<; (conjugate to m or 
to m'), which implies that they are symmetric under a Z(„ ^^-action and under a reflection, 
together giving a symmetry of D(„ ^). This means that if (n,£) > 1, then the configuration 
consists of k rings of (n,^) vortices in each. Typical deformed configurations with the 
correct symmetry in each mode for « = 3, . . . , 8 are shown in Figure 3.4. 
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(a) M = 4, £ = 2 



(b) n = 5,e = 2 




(c) n = 6J = 2 (d) (1 = 6, £ = 3 




(e) n = lj = 2 (f) ,1 = 7, f = 3 



Figure 3.4: Perturbations of the «-ring in mode £. These configurations are invariant under 
a subgroup isomorphic to D(„ ^j, and the grey lines in the centre of each represent the lines 
of reflection. The dotted figures are the regular n-gons. See section 3.4 for explanations. 
Continued on next page. 

In particular, when n is even and £ ~ n/2, the solutions have symmetry isomorphic to 
D«/2- Now in D„ sit 2 non-conjugate copies of D„/2^ one containing m, the other containing 

m', and since ^/"''^^ is fixed by m, the bifurcating solutions must have the symmetry T>„/2 
containing m. Consequently the bifurcating solution consists of a pair of regular «/2-gons, 
in general of different radii, staggered by Ijt/n as shown in Figures 3.4(a,d,k). 

When £ the D„-action factors through a Di.-action, and all bifurcating solutions 

have reflexive symmetry of m or m' (as seen from Figure 3.2). Now if k is odd, the resulting 
reflections are conjugate, so a configuration fixed by m will also be fixed by some conjugate 
of m'. This fact is illustrated in Figure 4(c) where n = 6, £ ^ 2, k = 3. Indeed, as D„ has 
n reflections while Dj^ has k, in the representation D„ D^, we must have {n,£) reflections 
in D„ that get identified, thereby fixing the same configurations. On the other hand, if k is 
even (as in n = 8, £ = 2, Figures 3.4(g,h), the nonconjugate m and m' in D„ have as images 
2 nonconjugate reflections in D^^, so the latter's fixed-point sets correspond to different 
configurations. 

Finally, whenever £ divides n, a perturbation in the -direction produces n/£ rings 
of .^-gons, in general of slightly different radii, and in the configurations with reflexive 
symmetry m the vortices in the different £-gons line up with the original «-gon. 
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(g) n = 8, f = 2 (Fix m) (h) « = 8, / = 2 (Fix m') 




(i) « = 8,^ = 3 (Fixm) (j) « = 8, £ = 3 (Fixm') (k)n = 8,£ = 4 

Figure 3.4: (continued) 



3.5 Degenerate critical points 

This section is dedicated to proving Theorem 6. We begin by presenting a criterion for a 
degenerate critical point to be a local minimum, then in sections 3.5.1 and 3.5.2 respectively 
apply the criterion to the cases where n is even and n is odd. 

Lemma 8. Let f be an analytic function defined on a neighbourhood of in M" with a 
degenerate critical point at 0, such that /(O) — 0. Write B — d^/(0), C = <3r'f{Q), D = 
dy(0). IfB is positive-semidefinite and if for all a G KerZ? \ {0}, b G M" we have 

Ca^=0, Da' + 6Ca\ + l>B\)^ >Q, 

then f is strictly positive on a punctured neighbourhood of 0. 

The symbol like Ca^b means 'evaluate the trilinear form C at a in 2 of its 3 arguments 
and at b in the 1 remaining argument'. In our application of this lemma, / will be the 
augmented Hamiltonian Hi (section 2.2), and for the point-vortex problem this is analytic. 

Proof. Suppose for a contradiction that ' (0) intersects every punctured neighbourhood 
ofO. 

Use the spUtting lemma (e.g. [19]) to write 

f{x,y) = Q{x)+h{y), 

where Q is a homogeneous quadratic form and /i is a function with vanishing 2-jet. Explic- 
itly, Q{x) — jX^Bx and y G KerB. Since y e KerB, the hypothesis Ca^ = implies that h in 
fact has vanishing 3-jet. Since the hypotheses of the lemma on the derivatives are intrinsic 
(invariant under change of coordinates), in the new coordinates x,y they become d^h — 
and (fh a'^>0 for all a e KerB \ {0}. 

Use the curve selection lemma [5, 9] to deduce the existence of an analytic curve 7(f) 
passing through along which / vanishes. In the new coordinates, write 7(f) = {x,y) = 
(^(f),Tj(f)); then {t)) + h{ri (t)) = 0. Expand Tj in Taylor series: 

^ (f ) = ^if + ^21'^' + ■■■, nit) = riit''+ J?2f + • • • , 
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where r is the order of the curve (at least one of £,i , Tji is nonzero). The leading terms of 
/o , Tj), all of which must vanish, are 

Qm)^\B^h''' + 3^1^21"-+' + ■■■ 

By inspection this forces = ■ ■ ■ — — 0- The coefficient of f'*'' is then B^^^^ + ^d^/i rjf, 
which must vanish. But the order of the curve being r and £,i — 0, we must have rji ^ 0, 

2 +i-c 



hence from the hypothesis Bi^,2_^i + ^d^^/i r\\ > 0, a contradiction. □ □ 



Consider now any Hamiltonian system on any symplectic manifold M, with Hamilto- 
nian H, symmetry group G, momentum map J : M — > g*, and suppose x e M lies on a 
relative equilibrium with finite stabilizer (possibly trivial). Define Tq := g - x n KerdJ. The 
symplectic slice M is then any G^-invariant complement to Tq in KerdJ. 

Proposition 9. If the hypotheses of Lemma 8 are satisfied with W replaced by Af and f by 
H, then the relative equilibrium is Lyapunov stable. 

Proof. It is enough to show that in the reduced space M^, jj. = J{x), the reduced Hamilto- 
nian admits a local extremum atx, cf. [17]. Since the action is locally free, J is a submersion 
near x. The tangent space of the submanifold J^' (ju) is To A/", and dj;^ {x) maps A/" iso- 
morphically to TxMj^, where x is the image of x in M^. Let / be a curve in through x. 
Then y lifts to a smooth curve in M tangent to M at x. The claim now follows because the 
hypotheses of Lemma 8 applied on are equivalent to the same hypotheses applied on 
TV. □ □ 

As the argument and calculations for the stability of the bifurcating points take distinct 
turns depending on the parity of n, we treat the even and odd cases separately. The even 
case is fairly easy, the odd case is much harder. 

3.5.1 n even 

The critical mode is ^ = n/2, we have i^t''/^) = a^'^l'^\ and both c and m act by multiplica- 
tion by — L This means that on this Fourier mode V„/2 = ^o:/"''^', Olg"^^''^, the augmented 

Hamiltonian Hi is an even function. We therefore expect, for generic families of functions, 
Hx restricted to V'„/2 to be equivalent to a family of the form 

fii(x,y) = ix'^ + Mx^+Z, 

where m is a parameter depending on ro, A. The + sign in front of the y^ term is jus- 
tified because the y-direction here corresponds to , whose eigenvalue from (2.6) is 
k:V/167i:>0. 

The D„-invariance of Hx helps us to figure out which terms arise in its Taylor series. 
For example, if / G V/, g € V^, and fg is invariant, then m — i. 

Lemma 10. For all a e V„/2 '^"'^ h^Af, we have Ca^b = 0. 

Proof. On the symplectic sUce expand Hx in Taylor series. Each term is invariant, in par- 
ticular the 3rd-order term Cx^ for x G J^. Given a, G Ve- (i — 1,2,3), the quantity Caia2a3 
lies in the tensor product V(^ V^j (g) Vf^, which contains invariant functions if and only if 
i'l ± ^2 ± ^3 = mod n. For ai = a2 = a G V„/2, Ca^b (necessarily invariant) can be nonzero 
only if b G Vq. But Vq n TV = {0}, implying that if b G TV and Ca^b is invariant, then b = 0, 
so that Ca^b = 0. □ □ 

It remains to calculate Da^ for a G V'„/2 ™ order to apply Proposition 9. The criterion 
for a bifurcation (3.6) reads, for even n, 

1+AM n^ 

(3.7) 



(1+Ar2)2 8(«-l)- 
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Put Zj = (ro + i-iyty-'^'j/" and f{t) =H^{zi,... ,z„). We shall expand 

/(') - - ji^ ^ £>»!!(>+ A - «.£ 

to the 4th order in f . 

Calculation At the bifurcation point, ro and X are related by (3.7) and CO ~ coq is given 
by (2.1). We are expanding 

/(O = L log {rl+t^- (4 -t'-)cos{2nk/n)) 

°^ KKn-l.A- odd 

"^"^^^ flog (ro + r)2 + log (ro-O^ 



32 TT 

-^fc^ (log(l +A(ro + r)2) +log(l +A(ro-f)')) 



2Vl+A(ro + f)2 l+A(ro-f) 

where • • • is a constant independent of ro, t, n, which will henceforth be ignored. Taking 
Taylor series in t of all of these terms to order 4 is simple, except for the first line, which 
comes out as 



n 



Stt ^ 327rr2 7687rr4 ^ ^ 

(up to an additive constant), thanks to identities akin to (2.4). The coefficient of in the 
Taylor series is then 



^^^(-(„-2)V+4(„-l)(l-Ar^f) 
which can be shown to vanish subject to the bifurcation relation (3.7). The coefficient of t 



-— ^ [(«-2)(«3 + 2«2-12« + 24)(t4 

+24(n - l)Ar,^(19A - 54AV;5 + 43Arg -4)] . 



Denote by T the term in square brackets. We wish to show that T > at the bifurcation 
point. Solving (3.7) for Arg yields 2 roots, substituting which into T in turn yields 2 
values, say Ti , T2 (functions of n). Write m = n — 2. A calculation (using Maple!) reveals 
that Ti + T2 is equal to 

128(ot+1) ,^ ^ 2560m + + A%32m^ + 5024m^ + I06\6m' 



[m^ - Am ~ 4) 

+ 15888m5 + 10778m'^ + 3266m^ + 177ot' 



while T\ T2 is equal to 

4096(m+l) , ^^^10 ^ ^ 7409ot*^ + 1044m^ + 39960m'^ + 855 12m^ 

(m2-4m-4)4 ^ 

+57332m4 + 25824m3 + 15 136m2 + 5376m + 576). 

In view of m ^ and n ^ 2, both are manifestly strictly positive, so that each of Ti and T2 
is indeed strictly positive. □ 
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Remark 11. Since the relative equilibrium is stable at the point of bifurcation, the result- 
ing pitchfork bifurcations are supercritical: the bifurcating relative equilibria are stable, 
and coexist with the unstable central one. Thus these stable bifurcating relative equilibria 
exist in a neighbourhood of the bifurcation point, Atq satisfying (2.11) with the inequality 
reversed. See also Figure 2.1. 

To persuade ourselves that this is a genuine pitchfork, we need to check the nondegen- 
eracy condition which is that the eigenvalues of the Hessian move through at nonzero 
speed with respect to the parameter Xr^ (or just A or ro separately). The expression (2.6) 
for the eigenvalues permits an easy check. 

As the mode that bifurcates is £ = n/2, the bifurcation occurs in the fixed-point space 
for the subgroup D„/2 as explained in section 3.4. The bifurcating solutions have D„/2" 
symmetry, i.e. consist of 2 regular n/2-gons at slightly different radii from the common 
centre, and these bifurcating solutions are stable, at least close to the bifurcation point. 



3.5.2 n odd 

There are two reasons why n odd is much harder than n even. First, KerB (degeneracy 
space) is 2-dimensional and its basis elements are less simple (namely ^''^^ rather than 
^(n/2)y ^QQQYid, the 3rd derivative contributions are nonzero and no analogue of Lemma 
10 holds. We proceed as far as we can with general odd n, and then specialize to numerical 
calculations for a few low values of n. 

The criterion for a bifurcation (3.6) reads, for odd n. 



1+l^r^ n + l 



(1+A 

The critical mode is £ = i (n — 1 ), and 



(3.8) 



K:=KerB=(a««-'V2\^(("-^)/2))cy(„_,)/2- 

We wish to apply Proposition 9, based on Lemma 8 with a G Vc and b G A/^. The calculations 
are simplified by the following observations. Recall the definition of Vj' from (2.8). 

Proposition 12. /. No cubic invariant exists on Vc, consequently Cs? = 0. 

Up to scalar multiple, there exists a unique quartic invariant on Vc, consequently 
Da"^ is a multiple of |a|^. 

Hi. Up to scalar multiple, there exists a unique cubic invariant of the form Ca^b with 
a£Vc and b G M, and invariance forces b G V[. 

Proof. Because j (n — 1 ) is coprime to n, the action of D„ on Vc is equivalent to the usual 
representation of D„ in the plane, though with an unusual choice of generator, cf. (3.2). 
D„-invariant functions on Vc are functions ofN = x^ +y^ and P = Re{x + iy)", cf. comment 
just before (3.5). Write a = {x,y). 

(i) As n ^ 5, this representation accommodates no cubic invariants, and as Ca^ must be 
invariant, it is 0. 

(ii) Likewise, the unique quartic invariant on Vc is A^^, so Da^ is a scalar multiple of 
^ |a|l 

(iii) If a G Vc and b G V„„ then Ca^b G V2 (n-i)/2+,« © V'2.(„-i)/2-mffi K,. For this to be 
invariant, we need m = or m— 1 . The former is ruled out by the assumption b G N, so 

bG7Vnyi=yi'. □ □ 

To understand better the invariant Ca^b, let x,y be as before on Vc and m,v be coor- 
dinates on V[, chosen so that m ■ = {x^—y) and m ■ (m,v) = (m,— v); in a nutshell 
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m ■ {z, w) — {z, w) in terms of complex variables z—x + iy and w = u + iv. We then have 
c ■ (z, w) = (c^"^ '^/^z, cw). The cubics of the form Ca^b are the real and imaginary parts of 
z^w, z?"w, |zpw. However, only the first of these is invariant under c, and only its real part 
is invariant under m. Thus, 

Ca^b = 7(z^w + z^vi>) ^ly{u(i?- -y^)-lvxy) 



for some value of 7 e K; explicitly 7=4 — ^ (Ca^b) — 7 

-.4 I /^„2k I ok2 - 



The key quantity Da + Ca b + Bb becomes, on completing the square. 



5|z|4 + 27Re(zM+i3|H^ = 5 



2 72 fi^_v2 



z + -^w 
o 



^kP. (3.9) 



Manifestly this is positive for all z, w ^ if and only if 5 > and j35 > 7^. The value of 
j3 is 

;2 



where as before a = 1 — Atq, cf. (2. 10). There remains the task of calculating 7 and 5. 

Calculation The awkward trigonometric expressions prevented us (and Maple) from 
reaching closed forms for 7 and 5. We therefore proceed to evaluate them numerically. 
In all these evalutations, ro is related to the parameter A by (3.8). For n = 7, of course 
A = and ro is arbitrary. 

„ 37.1 9.04 ^ 15.8 

= — , 7= 5 = —. 



''0 '0 '0 

j85 — 7^ = 504.7/rQ > hence the pentagon is stable. 

n = l : It transpires (Maple) that for n — lv^t have 

Computationally these numbers are correct to a high degree of precision — ^but we have no 
proof that they are rational multiples of 1/ Ttr^. At any rate it is certain that j35 — 7^ > 0, 
hence the heptagon is stable. This establishes Corollary 7. 

n = 9\ 

„ 6.9 13.8 ^ 182.4 

4 '•0 '•0 

Pd — y^ = 1075.60/rQ > hence the enneagon*^ is stable. 

«= 11: 

„ 12.0 29.7 ^ 555.1 

4 4 '"0 

j35 — 7^ = 5787. 6/rQ > hence the hendecagon'^ is stable. 
Conjecture 13. For all values of n, the n-gon is stable at the bifurcation point. 
Above, we have proved this for all even n and for n = 3,5,7,9,ll. 



'Nonagon' mixes Latin and Greek. 
' 'Undecagon' is anotlier Greco-Latin liybrid. 
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3.6 Bifurcations from the equator 

For all values of n > 3, the I ~ 1 mode 'bifurcates' at Xtq ~ 1, i.e. when the ring of 
vortices lies on the equator of the sphere. The momentum value at such a ring is fixed by 
all of S0(3): indeed, for the usual coadjoint equivariant momentum map for this problem, 
the momentum value is 0. Let us summarize from [14, Proposition 3.8] what bifurcations 
occur in this situation. On each near-zero momentum sphere we get, besides the regular 
ring of n vortices, the following configurations. 

n odd For each of the n planes through the poles of the sphere and containing one of the 
vortices, 2 configurations consisting of ^(n — 1) pairs and that 1 vortex on the plane; 
the vortices in each pair are each other's reflection in that plane. The notation in 
[14] is Ck{j{n — l)R, E), E referring to the single vortex on the plane, the R to the 
reflection pairs. 

n even In this case there are two distinct types of bifurcating solution, arising from the two 
distinct types of reflection in D„: 

i. For each of the planes through the poles and containing a pair of diametri- 
cally opposite vortices, 1 configuration consisting of \n—\ reflection pairs and 
those 2 vortices on the plane. The notation in [14] is Ch{{\n — 1)R, 2E). 

ii. For each of the planes through the poles and passing midway between adja- 
cent vortices, 1 configuration consisting of |« reflection pairs. The notation in 
[14] \^,Ch{\nR). 

4 What happens with other Hamiltonians 

Toward the end of sections 1.2 and 1.3, we met three options for families of Green's func- 
tions, which all agree for the plane A = 0. We have been opting for (1. 12). Here we sktech 
the conditions that guarantee the stability of the ring of identical vortices for the other two 
options (1.10), (1.11); the methods are the same as those of sections 2 and 3. 

4.1 Green's function G = log \z — wp 

Angular velocity The regular ring of n vortices with identical vorticities K rotates at 
angular velocity 

(0 = T-, 

Stt rl ' 

Unlike the expression (2. 1) this is not even in A. 

Stability The Hessian of the augmented Hamiltonian is the same as (2.2), except that A 
is changed to 

and the eigenvalues of the Hessian become 



while the expressions for Eg^' are as before. Also as before, the £ — [n/2j mode has the 
least eigenvalue, so the ring is stable provided Er"^^^ > 0. The criterion is 



(1+Ar2)2 2(«-l) 
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Compared with (2.11), for each n this new inequality is satisfied by a slightly narrower 
range of the effective parameter Xvq. The transition from stable to unstable still occurs at 
Atq of the same sign as for the previous Hamiltonian, and for A = the two Hamiltonians 
agree. Hence the Thomson heptagon is still stable. 



Bifurcations It seems likely that the bifurcations are of the same types as those explained 
in section 3; we have checked this for n = 5,7,9. 

4.2 Green's function G = log 

On the hyperbolic plane A < this reduces to the Hamiltonian adopted by Kimura [10]. 
For A = it is the standard Green's function on the plane, while for A > it corresponds 
to Green's function for the Laplacian on the sphere with 'counter- vortices'. The Hamilto- 
nian will model 2n vortices placed pairwise at antipodal points, each pair having opposite 
vorticities. Thus the 'ring' becomes 2 rings, one of n vortices of vorticity K near the North 
Pole, the other of n vortices of vorticity —K near the South Pole. In Laurent-Polz [12] these 
configurations are referred to as D„h{2R) when n is even and D„ci{R,R') when n is odd 
(in the former the rings are aligned, while in the latter they are staggered). The stability 
results of [12] are not directly applicable here, as he considers stability with respect to per- 
turbations of all 2n vortices, whereas we are considering a restricted class of perturbations: 
those preserving the antipodal pairing of the configurations. If a configuration is stable for 
Laurent-Polz, then a fortiori it will be stable for our setting. 

The calculations based on this option of Green's function get so cumbersome that the 
stability problem seems no longer tractable analytically. It seems likely that the results are 
similar to those in section 3, though the details of where the bifurcations occur will differ. 
We did calculate that the angular velocity ft) analogous to (2. 1) of the ring is 

^ \n - 1)(1 +Ar2) f 1 + i-Xry-') +2Xr^ ( 1 - (-Ar^)"" 



STZrl 1 - (-Ar2 

For A > (spheres) this can be deduced from [12, Proposition 3.11] 
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